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Mixed states of samples of spin s particles which are symmetric under permutations of the particles 
are described in terms of their total collective spin quantum numbers. We use this description to 
analyze the influence on spin squeezing due to imperfect initial state preparation. 
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I. INTRODUCTION 

Atomic samples for which external perturbations and 
mutual interactions between the atoms are completely 
symmetric under permutations of the atoms are conve- 
niently described in terms of collective variables. In par- 
ticular, if the atoms have two internal states relevant for 
the dynamical processes, each atom can be represented as 
a spin \ particle, and the full Hamiltonian of the system 
can be expressed in terms of the collective spin operator 
J = J2i A state vector of the system which reflects 
the permutation symmetry of n two-level particles can 
be expanded on the set of symmetric Dicke states | J° ) , 
where jo = n/2, and the quantum number m ranges from 
—jo, where all the atoms populate the states associated 
with a spin down single particle state, to j where all are 
in the spin up state. Intermediate values represent states 
containing all possible permutations of the atoms popu- 
lating either of the two spin i components, but with a 
definite total number, jo ± to, of atoms in the spin up 
(+) and down states (— ). 

Collective behavior of atomic ensembles is at the root 
of superradiance, and the Dicke states were introduced 
to account for atom-field interactions already many years 
ago. The past decade has provided numerous discussions 
and analyses of collective atomic dynamics in connection 
with atomic clocks and quantum information [Q, 0|, and 
recently a number of promising schemes has been pre- 
sented for tailoring of the collective atomic spin noise, 
relying on atom- light interaction j^, ||, collisional inter- 
actions in atomic samples [ 5 , 0. 0] , and on quantum non- 
demolition measurements [3,B|. 

Restricting calculations to the space spanned by the 
symmetric Dicke states implies a formidable simplifica- 
tion of the full 2™ dimensional problem to an n + 1 di- 
mensional one, and use of algebraic properties of the an- 
gular momentum operators further facilitate calculations 
of properties of the physical system. A large number 
of theoretical analyses therefore presents calculations re- 
stricted to this subspace. But, even if the assumption 
holds that the systems is invariant under exchange of 
atoms, not only the symmetric Dicke states will be pop- 
ulated: Imagine two two-level atoms, both described by 
the mixed state density matrix 



The combined state pi = P\® Pi of both atoms is also a 
mixed state. It is manifestly symmetric under exchange 
of the atoms, but if p ^ or 1 it is a mixture of all 
2 2 = 4 product states spanning the full Hilbert space of 
the pair of particles, and by a basis change we see that 
both the spin singlet and the three spin triplet eigenstates 
of the collective spin operator are populated. Unlike the 
spin singlet wave function, a density matrix component 
proportional to a projector on the spin singlet state is, 
indeed, invariant under exchange of the particles. 

Symmetric mixed states incoherently populate state 
vectors of different permutation symmetries. It is the 
purpose of the present paper to characterize this popu- 
lation of states of lower symmetry and to investigate the 
effect, e.g., on the achievements of spin squeezing and 
multi-particle entanglement of samples of atoms. For ex- 
periments, the situation with values of p close to zero 
or unity are particularly interesting as they represent the 
degree of polarization obtained, e.g., by optical pumping. 

In section |l| we derive some expectation values for the 
collective spin operators. In section III, we derive the 



probability distribution for states of different j and to, 
j < ioj f° r a collection of spin i particles in the state 
([j]). In section IV, we analyze the achievements of spin 
squeezing on such a collection of particles. In section [v| 
we discuss the case of s > |, i.e., atoms with more than 
two states. 



II. EXPECTATION VALUES OF THE 
COLLECTIVE SPIN OPERATORS 

Before determining the full density matrix of a system 
of uncorrelated atoms, it is instructive to see how much 
information can be derived from the fact that the atoms 
are uncorrelated. 

For a system of n uncorrelated atoms, we may de- 
termine the variances of the collective spin operators 
through 



AJ„ = nASi 



a 



x,y,z. 



(2) 



pi=pIT)(T| + (i-p)U)UI 



(i) 



In the case where the individual atoms are two-level 
atoms distributed according to (|l|), we obviously have 
(S 2 ) = | and (S 2 Z ) = \p + |(1 - p) = \, assuming 
h = 1. Since the distribution is invariant under rotations 
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around the z-axis, the x- and y-variances must be equal 
and wc find 



(3) 



Also, since the expectation value of S z is p — |, the vari- 
ance of J z is 



((S 2 )-(S 2 ) 2 ) = p(l-p)2j . 



(4) 



To find (J 2 ) we may again utilize the fact that the 
atoms are independent so that ( J 2 ) = 2i=y (^*) ' (^j) + 
Si (^i)' which in this case simplifies to 



(J 2 ) = (2p - l) 2 j (jo + 1) + 6p(l - p)jo- 



(5) 



In the completely polarized case, p = 1, the total spin 
equals jo = n/2 and (J 2 ) = jo (jo + 1)- Fot mixed en- 
sembles (J 2 ) is reduced, reflecting a distribution over 
states with different j, and in the completely unpolarized 
case, p = |, each atom contributes | to ( J 2 ) = n |. In 
the following section we shall determine the population 
of states with different j for any value of p. 



III. THE DENSITY MATRIX FOR A TENSOR 
PRODUCT OF MIXED STATES 

n spin i particles span a 2™ dimensional Hilbert space, 
7i, with basis vectors |toi, m„), m, = ±i. An alter- 
native representation is by means of eigenstates, | ^ ; A), 
of the collective spin operators J 2 (eigenvalue j(j + 1)) 



and J z (eigenvalue to), where A enumerates the degener- 
ate states within the eigenspaces, Ttj, m , of the J 2 and J z 
operators. 



A. Collective spin states for n spin | particles 

A complete description of the full set of basis states is 
given in Ref. ju|, but for completeness we shall derive 
our results by a formally simpler analysis. 

The | ^ ; A) states must span the full 2 n dimensional 
Hilbert space, and we start out by identifying the number 
of states with a given value of the quantum number to. 
Since to = J^. to^, the total number of atoms in the spin 
up state |t) is jo + m, and these states can be selected in 
a number of ways given by the combinatorial factor 



N m = 



n 

jo + m 



(jo +m)!(jo - m)\ 



(6) 



Due to the binomial formula (1 + 1)" = J2k (£)l*l (n_A,) 
we verify that J2 m N m = 2 n . 

States with given to values should now be grouped ac- 
cording to the values of j. Clearly j > to, hence we may 
begin with m = jo for which only a single state exists, 



and for which j must equal jo. The combinatorial fac- 
tor N m= j -i equals n (a single spin down particle can 



=30 

be chosen in n different ways) , and since only one of the 
states with m = jo — 1 can be obtained by the collec- 
tive lowering operator J_ acting on the state the 
remaining n — 1 states must have j = jo — 1. Applying 
the lowering operator J_ on all the identified m =jo — l 
states produces one state with j = jo and n — 1 states 
with j = jo — 1. The remaining N m —j —2 — n must there- 
fore have j = jo — 2. Continuing this argument provides 
the total number of states with j > m for any m > 0, 
and by symmetry we obtain the same number of j-states 
for negative m. In conclusion, the number of states with 
given values of j and m is independent of m and equal 



to D 



(1/2) 



N„ 



N„ 



=3+1- 



dim?i,-, m = D 



(1/2) 



2j + l 
jo+j + 1 



2jo 
jo + j 



B. The density matrix for a tensor product of 
mixed states 

Consider now the n-particlc density operator 



Pi 



(8) 



i.e., the tensor product of n identical two-state density 
matrices. Since any two-state density matrix may be 
diagonalizcd through rotations, we will only consider the 
diagonal form of p\ as given by ([j]) . 

As we intend to demonstrate the usefulness of the col- 
lective spin operator, we will only be interested in the re- 
duced density matrix, p™ d , obtained by taking the partial 
trace over the remaining degrees of freedom: 



red 
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A 



;A 



Pr. 



i -A 



(9) 



To determine the p-dependence of the matrix elements 
of p^ ed , we note that the vector | ^ ; A) may be expanded 
on the product basis, {|mi,TO2,... ,m„)}, so that (Q) 
only involves matrix elements of the form 



( 



m 1 ,m 2 , 



i I (Sm I 

, m n \Pf rn 1 ,m 2 , 



(10) 



where ^2 m 'i = m ' an d X) m » = m - Since p\ is diagonal, 
only matrix elements with m! i — mi are non- vanishing, 
and consequently p^ cd is diagonal in to. Writing the diag- 
onal elements of p\ as pi/ 2 + m s (\ _p)i/2-m s ^ ^ j g ev j(j en t 
that the value of any non-vanishing matrix element of 
the form (10) is p JO+ " l (l — p) JO_m . Since the expansions 
coefficients used to write p™ d in terms of these matrix 
elements are independent of p, we conclude that this is 
the full p-dependence of the matrix elements of p^ d . 

We now consider the case where p = |, so that the 
single particle density matrix, pi, is proportional to the 
identity. This implies that the full density matrix (||) is 
also proportional to the identity, or to be more specific 
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p n = 1/ dim7i in any basis. Consequently, since the en- 
tries of p^ ed are obtained by performing a partial trace 
as denned in (0), p™ d is diagonal, and the values of the 
diagonal elements are determined by the dimension of 
the relevant j, m subspaces, T~Cj, m . Denoting the diago- 
nal elements (populations) of p™ d byp^^, and using the 
dimension of Tlj jm found in ([?]), this gives us 



» 



dim7i 



j.»> 



dim7i 



= 2~ n D 



(1/2) 



for p 



1 



(11) 



Finally, for this form of p^ od at p = i to agree with the 
general p-dependence as found above, we find that p™ d 
must be diagonal in j for all values of p: Assume to the 
contrary that an off-diagonal element exists which is non- 
vanishing at some p-value. The required p-dependence 
implies that such a matrix element would also be non- 
vanishing at p = | in contradiction with (|TT|). Conse- 
quently such elements do not exist. 

In conclusion, 
plicit form of the diagonal elements is 



p r * d is completely diagonal, and the ex- 



(n) 
Pj.m 



2j + 1 / 2 jo 
io + j + 1 V Jo + j 



p> a+m (i-py 



(12) 



This surprisingly simple form of the full density ma- 
trix of an ensemble of identically distributed two-level 
systems immediately yields a host of interesting results. 

As an important example we may use the result to il- 
luminate the point that the states with less than full per- 
mutation symmetry arc populated. First, the marginal 
j-distribution, pj = 



, vfX^ is founc l t0 be 



Pj 



For p 



D 



(1/2) P 



2j+l 



(1-P) 



2j+l 



2p- 1 



-(p(i- p )y°- 



(13) 



Pj 



(2j + l)2- n D { ^ /2) according to dTTt) 



For convenience we shall assume p > | in the remainder 
of this section, and consequently it seems justifiable to 
approximate ( |13[ ) by dropping the term (1 — p) 2 i +1 in 
the numerator, thus obtaining 



Pj 



2 j + 1 p ( 2 j 
j +j + 12p-l Vjo+j 



Jo+i 



(l- p )j°-o. (14) 



We observe that the rightmost three terms constitute 
a normalized binomial distribution with an expectation 
value for j of j c = (2p — l)jo and a variance of a 2 = 
p(l — p)2jo. To justify the approximation made above in 
a more rigorous way, we may note that the omitted term 
corresponds to a similar binomial distribution centered 
at —j c - The approximation is thus expected to be valid 
as long as ~j c + a < 0, or equivalently p —\ > l/y/8 jo- 
The front-factor in the approximation ( [14] ) is slowly 
varying for j not too close to and it is close to unity at 
j = j c , which suggests that we may describe the marginal 
j distribution ( J13| ) by the binomial distribution alone. 
That this is indeed a good picture is illustrated in Fig. [I], 



where we also note that the limit of validity deducted 
above for the approximation (|i~i] ) seems to be correct, as 
the only peak where the approximation deviates visibly 
from the exact value corresponds to p — \ + 1 / y/8 jo ~ 
0.535 chosen at the validity limit. 

The variance in j may be approximated by the variance 
of the binomial distribution, a 2 , which according to (^) 
is equal to A J 2 . Since a 2 decreases with increasing p, we 
see that the j states close to the maximal value of jo are 
only populated at very high polarizations, in agreement 
with Fig. 0. In particular the total population of the 
symmetric Dicke states (j — jo) is given by 



1 



Pjo 



2p-l 



[P 



n+l 



-(1-P) 



(15) 



according to (|12|). To have pj g > Pj ~i we must have 
1 — p < 1/n, corresponding to an expectation of finding 
less than one atom in the spin down state. Note, however, 
that already at 1 — p < log2/n, we have half the total 
population in the symmetric Dicke states. 




FIG. 1: The distribution of the collective spin of 200 spin i 
particles at various polarizations. From left to right the peaks 
correspond to values of the polarization parameter p in the 
single particle density matrix (|l|) of 0.535, 0.65, 0.75, 0.85, 
0.95 and 0.99 respectively. The solid line is the exact value 
of pj according to (113), the dashed line corresponds to the 
approximation (^) which cannot not be discerned from the 
exact values at large p- values. Note that the p — 0.99 line 
ending in mid-air is not a technical artifact, but illustrates a 
significant population of the symmetric Dicke states as given 
by (§). 



IV. THE EFFECT OF SPIN SQUEEZING ON A 
COLLECTION OF SPIN | PARTICLES 

The variance of the spin components orthogonal to the 
mean spin of the state (JsJ) obeys AJ X AJ y — jo/2, which 
exceeds the Heisenberg uncertainty limit, | \{[Jx, Jy]}\) 
by a factor of 1 2p — 1 1 1 , which approaches unity in the 
highly polarized cases (p = 0, 1) as expected. 

It is interesting to study the effect of spin squeezing op- 
erations on such states, since they might serve as models 
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of imperfectly prepared initial states on which squeezing 
operations are performed. Note that obtaining a spin 
squeezed state indicates that the initially uncorrelated 
atoms described by (||) have become entangled JllJ . 

Spin squeezing may be performed in a number of ways, 
and we will discuss the two unitary squeezing methods 
analyzed in Ref. jl2| in the case of initially pure distribu- 
tions. Both of these approaches make use of a Hamilto- 
nian which is quadratic in collective spin variables, i.e., it 
involves interactions among the particles. The descrip- 
tion of the evolution of the states would be immensely 
complicated in the product state representation of the 
many spins where, in fact, very entangled states are pro- 
duced. In the collective spin representation, however, the 
action only has to be computed for each total spin com- 
ponent, which is a much easier task. 



A. One-axis twisting 

A one-axis twist about the y-axis, as mediated by the 
Hamiltonian \ Jy > will squeeze the variances of the J x and 
J y spin components of an initially maximally polarized 
state as illustrated on Fig. ^ (the dotted ellipse): The 
variance of the spin along one direction in the xy-plane 
is increased, while the variance along the perpendicu- 
lar direction is decreased to a value below the standard 
quantum limit jo/2. 

As we have tried to illustrate by an example on the 
same figure, the same is qualitatively true in the case 
where the initial state is not a coherent spin state, but a 
"sufficiently polarized" mixed state. 

Wave function components of definite j show a reduced 
variance 



A J? = I I 1 - 



2j-l 



A] + B]-A 3 



(16) 



of the spin component Jg 
i arctan-? 1 , As = 1 — cos 



cos 9 J v — sin 9 J x , where 9 



2j-2 



/t. 



B* 



= 4sinf cos 2 ^ 2 § 



2 JT) 

and (i = 2\t, fOq |. We note that both the optimum in- 
teraction time and the direction 9 depend on j. For a 
distribution of j-values this will deteriorate the squeez- 
ing. Furthermore, the results in ]T^ | assume an initial 
state with m = j = jo, whereas we start out with a 
distribution of m- values. 

Assuming that p is so close to unity that we may use j c 
and a from section [II B to estimate the center and width 



of the j-distribution, the relative variation of the optimal 
fi over the width of the j-distribution is approximately 
A/x /// ~ ~p)/jo- Consequently, for j > 1 it is 

possible to choose a value of \i that is nearly optimal for 
all populated j subspaces. According to ( |l2| ) and (13), 
we find that within each j-subspace, the fraction of the 
population in the uppermost m level, p^ , is larger than 

Similarly the direction of the squeeze axis does not 
seem to vary significantly over the populated j subspaces, 



as is also supported by Fig. || showing results for 15 atoms 
with p = 0.9. Looking at the distributions within the dif- 
ferent j-subspaces (dashed ellipses) , one should note that 
the outermost ellipse corresponds not to the maximal j 
but to j = jo — 1, which has the greatest population, jo 
correspond to the second-largest ellipse and we can see 
how the squeezing axis turns clockwise with increasing 
j. When we perform the sum over all j-values we obtain 
the result illustrated by the solid ellipse in Fig. ||. The 
state is significantly squeezed, but less than in the ideal 
case indicated by the dotted ellipse. As anticipated, the 
population of the states with non-maximal m within each 
subspace does not appear to have serious adverse effect 
on the squeezing process. 





— 
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FIG. 2: The distribution of the spin components orthogonal 
to the mean spin after performing a one-axis twisting squeeze 
on a system of 15 spin i particles with initial single particle 
density matrix given by (hi) with p = 0.9. The squared length 
of the minor and major axes of the ellipse drawn with a solid 
line equal the variance of the spin-component in these direc- 
tions. Additionally the direction of the minor axis correspond 
to the direction of minimal variance. The dotted ellipse shows, 
for comparison, the result of applying an optimal squeezing 
to an initially pure system (p — 1), and the dashed ellipses 
illustrate the contributions from the various j subspaces in 
the mixed case as discussed in the text. 



B. Two-axis countertwisting 

We analyze the time evolution under the two-axis 
countertwisting Hamiltonian, H c t = xj2i {j\ — J 2 .) . 
One advantage of the two-axis countertwisting Hamilto- 
nian is that the squeezing is always in the same direction. 
A disadvantage is that it cannot be solved analytically, 
and we will proceed according to the Holstcin-Primakov 
approximation which is only valid for highly polarized 
samples and not too strong squeezing. Consequently we 
shall restrict the validity of our results to lowest order in 
1 — p and assume jo to be much larger than unity. Us- 
ing results from Jl3[ we find the variance of the squeezed 
component, J y , of J after squeezing to be 



,m~(j(i + l)-m 2 )e- Ai , (17) 
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where A = 4%i. In the fully polarized limit (p = 1) 
we have AJ y 2 = exp(—Xjo) jo/2, and consequently Sq — 
exp(Aj'o/2) is the pure state squeezing factor. 

It is a straightforward task to perform the sum over 
m in ([l?]). Doing this we obtain two terms inside the 
j-sum, one of which may safely be dropped under the 
assumptions already made (p > | and jo 3> 1). Keeping 
the remaining term we find 



AJ, 



2 ~ e -Ajo [ £ 



2 jo 



E — 



Pa / (2p-l) 2 
2 ' 2io 



3 + Jo \jo + j 



(18) 



where we have used the asymptotic approximation, jo ^> 
1, extensively and introduced p\ a,sp\ = p/(l + (exp(A) — 
1)(1 — p)). The three last terms of the second line are 
recognized as a normalized, rather narrow, binomial dis- 
tribution. Approximating this distribution with a 5- 
function centered at j = (2p\ — l)jo we find 



v 2 \pj \2p-l 

To first order in (1 — p) this is equal to 



(19) 



AJ U 



e-^Jo ( 1 + (e A_ 1)(2io _ 3)(1 _ p) ) 



Based on this result we may determine the effect of 
the non-purity of the initial distribution on £r: the ratio 
between the frequency uncertainties obtainable by Ram- 
sey spectroscopy on squeezed and unsqueezed samples, 
as introduced by Wineland et.al. [Q. According to this 
reference we have 



AVI(J. 



(Ac) (AJ y /\(J z )\) ' 



(20) 



where the terms in the denominator refer to the un- 
squeezed sample. Since (AJ 2 ) Q = jo/2 is independent 
of p, and (J z ) = ((J z ))q to first order in A, we find that 
to lowest order in (1 — p) and A we have 

&«-(l + 2(l-p)logao), (21) 
so 

that is: to first order in 1 — p we get a correction to 
the quality of the squeezed ensemble as compared to the 
case of performing the same squeezing operation on an 
initially perfectly polarized ensemble. 

Note that we have only examined how imperfect initial 
state preparation affects the Ramsey spectroscopy; for 
an analysis of the significance of decoherence during the 
measurement see Refs. 14, p"5|. 



V. ATOMS WITH MORE THAN TWO LEVELS 

One of the virtues of the combinatorial analysis of col- 
lective spin distributions is that it may be generalized to 
the case where the atoms have spin s > 4. 



A. Collective spin states for n spin s particles 



Proceeding as in section III A we find that the dimen- 
sion of the | -j, ) subspace is 



D 



W _ 



jo + j 



2 s 



n 

io+j + i 



(22) 



2 s 



where the generalized binomial Q^J . denotes the number 
of ways of placing m indistinguishable balls in n distin- 
guishable urns so that no urn contains more than k balls. 
Note that this definition agrees with that of the ordinary 
binomial for k = 1. 

The main problem with this approach is that the gen- 
eralized binomials are not as analytically accessible as the 
binomials. Their numerical values may, however, easily 
be obtained through the recursion relation 



E 



71-1 

m — i 



(23) 



and their analytical properties are to some extent acces- 
sible from the generating function: 



E 



m 



1 + t + . . . + t* 
l _ t k+i- 



1 -t 



(24) 



B. The density matrix for a tensor product of 
mixed states 

In order to represent the tensor product p n = pf n in 
the basis of eigenstates of the collective angular momen- 
tum, we have to develop this basis and to identify the 
matrix elements of p n . It turns out to be advantageous 
to perform these two tasks in parallel by sequential cou- 
pling of angular momenta of 1, 2, . . . , k spin s particles, 
and determination of the appropriate density matrices 

Pl,p2, ■■■ ,Pk- 

Let us introduce the collective spin of the first k par- 
ticles as 



(25) 



We denote by {\ji,j%, ■ ■ ■ ,jk, mk)} the simultaneous 
eigenstates of (J*- 1 -*) 2 , . . . , (J^-*) 2 and with eigen- 
values ji(ji + 1),... ,jk{jk + 1) and mfe respectively. 
These states are obtained by sequential coupling of the 
atomic spins: by coupling the first, j\ = s, and second 
atom we obtain a total spin j2, which is subsequently 
coupled with the third atom to the resulting spin j'3, etc. 

The density matrix p\ is our known input. We assume 
that m s labels the basis vectors in which p\ is diagonal 
with eigenvalues p m „ (note that m s does not need to be 
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associated with any real angular momentum or spin of 
the physical system, the spin representation is convenient 
to represent any (2s + l)-level atomic system). 

If we assume that the density matrix, p k , of k atoms is 
known in the basis of angular momentum coupled states, 
{\jlihi ■■■ ijk,Tnk)}, we can recursively construct the 
density matrix of k + 1 atoms: 



Pk+i = ^2 \a' k ,m k ) (a k ,m k \ 

a k ,a' k ,m k ,m s 

(a' k ,m k \p k \a k ,m k ) p„ ls , 



(26) 



where we have utilized that p k is diagonal in m k as 
argued in section III B, and we have introduced a k 



as a shorthand for the intermediate angular momenta 
ill 32 1 ■ ■ ■ ijk- It is now just a matter of angular momen- 
tum coupling to obtain from ( p6] ) the density matrix pk+i 
in the basis of total angular momentum states of the k+1 
particles, specifying j k +i, Wfc+i an d all the intermediate 
angular momenta a k . 

The bookkeeping seems to be a demanding problem, 
but as in Eq. (|) we are only interested in the reduced 
density matrix retaining the properties of the total spin 
of the k+1 particles. Hence we trace over a k , i.e., we 
add incoherently the contributions from subspaces with 

(k) 

different intermediate jy : Let p\ ' denote the diagonal 
elements of the reduced density matrix p k ed , that is, the 
probability that k atoms have a total angular momentum 
of jk and projection m k . We can then write the matrix 
elements of the reduced density matrix for k + 1 atoms: 



the combined state of these particles is 



2 I 2 \ 

3 I 0/ 



1 |0\ 
3 I 0/ i 



(28) 



which clearly yields a non-diagonal reduced density ma- 
trix. 

Note that according to (12) we actually have the 
stronger selection rule Aj = < 2s for s = \. 

In closing, we present a simple case of a mixed state 
of a system of uncorrelated atoms with more than two 
levels, which might serve as an approximation of e.g., an 
optically pumped system. We assume the single atom 
density matrix, p\, to be diagonal, and of the form 



Pn 



= L e -(3m s 



(29) 



where Z is a normalization constant depending on the 
parameter (3. By induction through (p7|), or by the argu- 
ment that we used to deduce the p dependence in section 
IIIB, we find that with this form of p\, the parameter 
dependence of the reduced density matrix, p^ ed , must be 
Z~ n exp(— /3m). By the same arguments as in section 



IIIB we conclude that the reduced density matrix is di- 



agonal and 



where Dj is given by (|2 



„(") _ n( s ) 7 -n -f3m 



(30) 



3k+i 
n'k+i 

E 

jk,m k ,m 



red 
Pk+1 



Jk + l \ _ 



(27) 



m 'k+i 



Jk s 
rrik m s 



Jk s 
rrik m s 



jk + l \ „ 

"ifc+i /Pj k ,m k Prn s , 



where the brackets, (ii]° k m s \ m k +\)i denote Clebsch- 
Gordan coefficients. An interesting consequence of 
the relation (|27j ) is that for k less than n it is 
sufficient to determine recursively the diagonal ele- 
ments of the reduced density matrix, p T k ed : 



E 



jfc_i,mfc_i,m 3 



K 



mfc_i m s 





m k 



k ) 



2 (fe 

p 



j k - 1 ,m k 



.! Pn 



J k ,m k 

In the 



last step we obtain p™ d which also has coherences be- 
tween different collective spin states. According to the 
selection rules of angular momentum coupling this im- 
plies that the reduced density matrix must be diagonal 
in m n and band diagonal in j n with values ranging over 
\Aj n \<2s. 

As an example, consider two spin 1 particles which are 
both in the pure state | J ) . In the collective spin basis 



VI. SUMMARY 

In conclusion, we have reiterated the convenience of 
collective spin operators for analyzing the behavior of a 
system of atoms. 

We have established the explicit form of the reduced 
density matrices of ensembles of two and multi-level 
atoms with a given single particle distribution, and we 
have used these explicit distributions to predict the ef- 
fect of non-perfect sample preparation on the effective- 
ness of spin squeezing by one-axis twisting and two-axis 
countertwisting. 
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